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Abstract. Semiconductor-superconductor hybrid systems are promising candi- 
dates for the realization Majorana fermions and topological order, i.e. topologi- 
cally protected degeneracies, in solid state devices. We show that the topologi- 
cal order is mirrored in the excitation spectra and can be observed in nonlinear 
Coulomb blockade transport through a ring-shaped nanowire. Especially, the exci- 
tation spectrum is almost independent of magnetic flux in the topologically trivial 
phase but acquires a characteristic h/e magnetic flux periodicity in the nontrivial 
phase. The transition between the trivial and nontrivial phase is reflected in the 
closing and reopening of an excitation gap. We show that the signatures of topo- 
logical order are robust against details of the geometry, electrostatic disorder, and 
the existence of additional subbands and only rely on the topology of the nanowire 
and the existence of a superconducting gap. Finally, we show that the coherence 
length in the nontrivial phase is much longer than in the trivial phase. This opens 
the possibility to coat the nanowire with superconducting nanograins and thereby 
significantly reduce the current due to cotunneling of Cooper pairs and to enhance 
the Coulomb charging energy without destroying the superconducting gap. 



PACS numbers: 74.25.F-, 85.35.Gv, 74.78. Na, 74.20.Rp 
1. Introduction 

Topological phases are quantum phases which cannot be described by a local order 
parameter. Instead, the defining characteristic of topological phases is a pattern 
of long-range quantum entanglement which is called topological order [TJ [21 [3J. 
One characteristic property of topological order is the dependence of the ground- 
state degeneracy on the topology of the manifold on which the system is defined 
[2 21 IH). This degeneracy on manifolds might also serve as a starting point for a 
general classification of topological phases of strongly correlated quantum matter, 
complementary to the topological band theory which is based on single-particle states 
and cannot be easily generalized to correlated systems [6j [7] . Recently, there is much 
interest in topological phases [8, 9 due to their possible application in topological 
quantum computation. However, these phases are also of fundamental scientific 
interest for their ability to support exotic quasiparticle (QP) excitations with abelian 
and even nonabelian quantum statistics. 
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One particularly interesting class of topological phases are topological 
superconductors, which have been predicted to host Majorana bound states [^I HUIUTI 

[IimmimmHTlHillHlHD]. ThePx+iPv su P erconductor ( SC ) for s P inless fermions 
is a prototype system for topological SCs. Depending on the chemical potential, the 

ground state of the p x + ip y SC is realized by the weak or the strong pairing phase, 

which can be distinguished topologically. In particular, the grand canonical ground 

state of the weak pairing phase on the torus depends on boundary conditions (BCs) for 

each of the two primitive directions [H [5J . Here, the ground state with only periodic 

BCs is special and shows an odd parity, while the three ground states with at least one 

antipcriodic BC are characterized by an even parity. In contrast, the strong pairing 

phase and also the ordinary s-wave SC on the torus possess a fourfold degenerate even 

parity ground state [5J. 

In this paper, we consider a quasi one-dimensional ring shaped SC nanowire 
and demonstrate that essential aspects of the above described topological degeneracy 
on the torus carry over to this simpler geometry. We focus on a regime in which 
the quasiparticle gap A is larger than the single-particle level spacing d. In the 
Coulomb blockade regime, the mean particle number N and the parity P are fixed 
by the charging energy E c > A and the degeneracy of grand canonical ground states 
is reflected in the excitation energies, which can be observed in nonlinear Coulomb 
blockade transport [2I1I22]- The lowest excited state above the ground state of a trivial 
SC with even parity involves two Bogoliubov QPs and thus breaks a Cooper pair, 
incurring an excitation energy 8E w 2 A, which is essentially independent of magnetic 
flux [23]. In contrast, the ground state for odd parity always has one Bogoliubov QP, 
and hence the lowest excited state involves both annihilating and creating a Bogoliubov 
QP which costs the excitation energy 8E d 2 /A <C 2A. For nontrivial topological 
SCs the situation is very different. Here, ground states without unpaired particles at 
the Fermi energy have odd parity for periodic BC, and even parity for anti-periodic 
BC. Therefore, the excitation energy SE oscillates between d 2 / A and 2 A as function 
of magnetic flux with period h/e which is doubled as compared to the case of a trivial 
SC [53] • This connection between the ground-state degeneracy on manifolds with 
nonzero genus and the h/e flux periodicity of ring structures demonstrates that these 
properties are a general consequence of topological order and that nonlinear Coulomb 
blockade transport is a suitable tool to investigate topological order. 

Recent experiments show evidence for Majorana fermions and topological 
superconductivity [231 123 [23 [371 I2E] in semiconductor (SM) nanowires with strong 
Rashba spin-orbit coupling in a magnetic field and proximity coupled to an s-wave 
SC [2S1 EH EH E2]- Despite the effort which has been invested, clear experimental 
signatures of the unconventional nature of the superconducting state are still missing. 
Therefore, additional detection schemes have been suggested such as the periodicity 
of the Josephson effect [23 [Ml EH [35] > tunneling spectroscopy ESI EH ESI EH 
HO] , interferometry [HI 22], transport experiments [33J 231 SSJ, and coincidence 
measurements [46 j . We here propose another experiment which directly investigates 
consequences of topological order on a nontrivial manifold. For this purpose, we use 
the Coulomb energy as an instrument to prescribe the parity of the hybrid system and 
thus to observe the above discussed ground-state degeneracy. Our analysis is based on 
the identification of the pfaffian Z2 invariant Q for Hamiltonians in class D [33] with the 
parity of the grand canonical ground state. Thus, we use this key piece of information 
about the grand canonical ground state to construct two classes of states with parity 
Q and —Q, where the class of states with parity Q (— Q) contains all eigenstates with 
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an even (odd) number of QP excitations. We find two types of excitation spectra 
which display trivial or nontrivial superconductivity depending on parameter values. 
The transition between the different topological phases is characterized by the closing 
and reopening of an excitation gap. As these findings only rely on the existence of a 
superconducting gap A > d and the S 1 topology of the system, the excitation spectra 
are robust against disorder, spatial variations of the superconducting pairing potential, 
geometry details, and the existence of additional transverse subbands. 

This paper is organized as follows: In section [2j we introduce the model system 
and the proposed experimental setup. We continue in section [3] with a brief review of 
the results in Ref. [23] for single-band SM nanowires and study the robustness of these 
results against details of the geometry, electrostatic disorder, and local variations of 
the superconducting order parameter. In section [4] we make a departure from the 
case of strictly one-dimensional nanowires and consider the experimentally realistic 
case of quasi one-dimensional nanowires. In section [5j we compare the current for the 
single-electron tunneling with the current due to cotunneling of Cooper pairs which is 
the most relevant transport channel competing with sequential tunneling of electrons. 
We summarize our results in section |6] 

2. Model system 

We consider a quasi one-dimensional SM nanowire with strong spin-orbit coupling 
which is proximity coupled to an s-wave SC. The nanowire forms an annulus in the x- 
y plane with radius R and radial extension L± -C R. We assume a strong confinement 
in z-direction, i.e. the extension perpendicular to the plane of the annulus L z <C Lj_, 
such that only the lowest subband with momentum in z-direction is occupied. This 
hybrid system is separated from a back-gate by a thin insulating layer and weakly 
tunnel-coupled to source and drain electrodes with potentials ±eV/2. Assuming a 
strong capacitive coupling between the nanowire and the SC, the total number of 
electrons in this system is determined by the Coulomb Hamiltonian 

H c = E C (N W + N sc ) 2 ~ eV G (N w + N sc ), (1) 

where E c denotes the charging energy and Nw (Nsc) the number of excess electrons 
in the SM (SC) attracted by the back-gate. Varying the gate potential eVc allows to 
change the total electron number N — N\y + Nsc m discrete units. A current through 
this island involves changing the electron number from N to N ± 1 and creating or 
annihilating a QP excitation. Thus, resonances in the differential conductivity appear 
when the condition 

eV/2 = E N±1 - E% (2) 

is satisfied, where Em is the total energy of a state with N electrons and 
the respective ground-state energy. The spacing between the resonance peaks is 
independent of the charging energy E c and displays the excitation spectrum for fixed 
particle number, 

SE N = E N - E%. (3) 

In our analysis, we assume that both the charging energy E c and the QP gap A$c 
in the SC are larger than the effective gap A e g in the SM. Hence, for small voltages 
eV < E c ,Asc all electrons in the SC are paired and unpaired electrons as well as 
breaking of Cooper pairs can only show up in the SM. 
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We describe the low-energy physics of the nanowire by the lattice Hamiltonian 

h = h sm + h sc [H M El 132] with 

H S U = { ~ trr ' + ^ r ' r ' ( Ez °™ - V + V r ) }4 CT c r'o- 

r,cr 

~ 2a XI (4<T^aff c r+««S - 4*Va& c r-8 yS ) , (4) 

where the operator c\ a (c rCT ) creates (annihilates) an electron at site r with spin a 
and mass m. The first term describes hopping on a simple square lattice with lattice 
parameter a, i r . r +<5 = to = ti 2 /2m 2 a 2 for the nearest-neighbor lattice vectors d, and 
i r r = —2^o- The second term in Eq. Q contains the chemical potential fi, the 
electrostatic disorder potential V r , and the Zeeman energy splitting Ez = gfisB/2 
due to the magnetic field in z-direction. The last terms in Eq. Q represent the Rashba 
spin-orbit coupling with spin-orbit velocity a, and a s are the Pauli spin matrices with 
s = x, y, z. By coupling electrons with opposite spins, the spin-orbit coupling creates 
two helical bands with spin rotating in the x-y-plane. The magnetic field tilts the 
spin direction out of the x-y plane and removes a level crossing at zero momentum by 
opening a spin gap. 

The proximity coupling between the s-wave SC and the nanowire is described by 
the effective s-wave pairing Hamiltonian 

H SC = { A r C U C U + A ^riCr T ) (5) 

r 

with pairing potential A r . The pairing potential modifies the dispersion relation in a 
crucial way by inducing two excitation energies. It opens an effective pairing gap at 
the Fermi surface, and reduces the spin gap at zero momentum. Depending on Ez, 
fx, and A, the SM nanowire shows two topologically distinct superconducting phases 
which are separated from each other by a topological phase transition [33 131] ■ 
For one partially occupied subband and E z < A 2 + fj, 2 , the nanowire is in the 
topologically trivial phase. By tuning the chemical potential via the gate voltage 
or the Zeeman energy via the magnetic field, the nanowire can be tuned across the 
phase transition, which shows up as the closing and reopening of the QP excitation 
gap. In the topologically nontrivial phase, which is reached for E 2 Z > A 2 + [i 2 , such 
a nanowire supports a pair of zero-energy Majorana bound states located at the ends 
of the wire. 

We here consider a closed nanowire without end points and thus without 
Majorana bound states. However, the unconventional nature of the topologically 
nontrivial superconducting phase shows up in the doubling of the magnetic flux period 
of the excitation spectrum from h/2e to $o = h/e [23 . In our analysis, we incorporate 
the magnetic flux $ through the azimuthal vector potential A = <fre v /2irR and the 
Peierls substitution with 



if S^ s A(r')dr' 



A -> Ae iq ' r , 



(6a) 

(66) 
(6c) 
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where q is the Cooper pair wavenumber. We determine the Cooper pair wavenumber 
by minimizing the ground-state energy of the nanowire with respect to arbitrary 
vectors q. This is equivalent to determining q by minimizing the Ginzburg-Landau 
free energy for the s-wave SC coupled to the nanowire 23J . For the strictly one- 
dimensional ring-shaped nanowire with azimuthal q = qe v this demands that q is the 
integer nearest to — 2$/$o- 

We diagonalize the Hamiltonian by defining the Bogoliubov QP operators oli — 
Ylr a ( u raic rC r + v r <yic\a) i wnere {1} 1S a complete set of QP quantum numbers. This 
yields 



H=^E l a\a l +E GG (7) 
i 



with Ei > and the ground-state energy 

lj2( 2t 0-V + Vr)- (8) 



e gc = -\Y, E i ' •> 

I 



The corresponding ground-state electron number is given by the expectation value of 
the particle number operator N — ^2 rcr c \ a c ra hi the state where all QP levels are 
empty. Rewriting the particle number operator in terms of QP operators and taking 
the expectation value with respect to the ground state, we find 

iVcc^El^l 2 - ( 9 ) 

I-CT I 

The parity of the grand canonical ground state is determined by the pfaffian Z 2 
topological number 

Q = « , (10) 
Vldet (Wt x ) | 

where T~L denotes the Bogoliubov-de Gennes Hamilton matrix in the basis 
(cL ., cL , c r f, c r i) and t x denotes the Pauli matrix acting on the particle-hole space 
(331 147]. Here, the topological number Q = +1(— 1) corresponds to the even (odd) 
parity of the grand canonical ground state. From the grand canonical ground state 
with parity Q we construct two classes of states with parity P by creating N qp QP 
excitations. The parity P of these states is determined by 

P=Q.{-\) N *», (11) 

i.e. depending on whether iV qp is even (odd), P = Q (P = — Q)). 

Since the large Coulomb energy enforces a well-defined parity and mean particle 
number in the nanowire, we Legendre transform both classes of states into a pseudo- 
canonical ensemble with N qp QPs, mean particle number N, and parity P. This 
ensemble contains states with energy 

E{{l i },N,P]=E GG + Y / Ei j +liN, (12) 
i=i 

where the chemical potential fi is determined by the constraint 

N = N GC (ri (k^(M)l 2 - k^(M)| 2 )- (13) 

j = l re 



Robustness of Topological Order in Semiconductor-Superconductor Nanowires 



G 



Here, {U : i — 1, ...,N qp } denotes the set of quasiparticle excitations. For 
noninteracting systems, this method is fully equivalent to the wave function based 
technique used in Ref . [33] . There, the authors defined two classes of grand canonical 
ansatz wave functions with even and odd parity as tensor products of generalized 
BCS wave functions [55]. The ground state was then determined by the unbiased 
minimization of the energy expectation value and the lowest excited states where given 
by pairwise creation of Bogoliubov QPs. However, while the wave function method is 
defined on the many-body Hilbert space, the technique used here only relies on the 
pfaffian Z2 invariant, the QP energies, and the corresponding QP wave functions and 
is thus more suitable for large systems without additional symmetries. 

In recent experiments, it has been shown that both InSb and InAs are suitable 
SM materials due to a large g-factor and strong spin-orbit coupling '24, 25, 26 ] |27 | 128] . 
In the experimental situation, the confinement energy in transverse direction is the 
largest energy scale so that only a few subbands are partially occupied. It is useful 
to express the kinetic energy and the spin-orbit coupling in terms of a characteristic 
energy e so = ma 2 /h 2 and a spin-orbit length l so = fi? /ma. Typical values for these 
materials are e so — 0.1 meV, Z so — 100 nm, and <7/x^/2 — lmeV/T. Thus, we 
find with R = 0.5 pm, A = 0.5 meV, and E z — 1 meV, an effective pairing gap of 
A ef j ss 0.2 meV and the level spacing at the Fermi energy of d = 0.08 meV. To ensure 
sequential single-electron tunneling through the hybrid system, we consider the case 
E c ^> A e ff. Here, a relative charging energy between the SM and the SC of the order 
of 1 meV together with a pairing gap A$c = 2 meV in the SC [24] would reduce A e g 
by 20 % [48] and is thus neglected. In the following, we vary the magnetic field B 
in discrete steps with the magnetic flux always being fixed modulo $01 such that the 
only effect of B is a change of the Zeeman energy Ez ■ 



3. Single-band Hamiltonian 

We begin our analysis with a perfectly one-dimensional nanowire of width Lj_ — > 
which has been studied recently in Ref. [23J. There, the authors considered the 
rotational symmetric continuum version of the model introduced above, i.e. spatially 
constant superconducting pairing A r = A and vanishing electrostatic disorder V r = 0. 
Before discussing the relevance of geometry details, electrostatic disorder, and spatial 
dependence of the pairing potential, we recapitulate the main results of Ref. [23]. To 
illustrate the main results in that Letter, we will here focus on the nondisordered case 
(7 = 0, black line) in figure [2] The disordered case will be treated in section 3.2 The 
plot shows the lowest excitation energy En — Ef^ as a function of Zeeman energy with 
N chosen such that /x w 0. As discussed above, at E 2 Z = A 2 + p? the Zeeman energy 
drives the nanowire through a topological phase transition with the trivial phase for 
Ez < A and the nontrivial phase for Ez > A [T5] . 

Our findings for Ez < A are characteristic for s-wave superconductivity in 
metallic nanograins |21[ I22j . For even parity, the excitation spectrum shows a 
superconducting gap 2A e g since all excited states contain two Bogoliubov QPs which 
corresponds two breaking one Cooper pair. In contrast, the ground state for odd parity 
always has one Bogoliubov QP and therefore the spectrum is qualitatively independent 
of magnetic flux and determined by the single-particle level spacing as 

6E = E(n = i)-E{n = 0)n^- (14) 
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for QP energies E(n) — \Jn 2 d 2 + A 2 S where n counts the energy levels relative to the 
Fermi energy. Similarly, a variation of magnetic flux by $o/2 changes E(n) on the 
order of c? 2 /A c g. 

In the topologically nontrivial phase for Ez > A, the parity effect is very different. 
Here, the excitation energies depend on both electron parity and magnetic flux. In 
Figs, ga) and (d) we find a QP excitation gap 2A G ff since two Bogoliubov QP 
excitations are required and thus a Cooper pair needs to be broken. In contrast, 
the excitation energies in Figs. |2jb) and[2jc) reflect the single-particle level spacing 
as <i 2 /2A e ff since always one unpaired particle is located near the Fermi surface. As 
shown in Fig. [2^d) , the characteristic signature of the topological phase transition at 
Ez ~ A is the closing and reopening of the QP excitation gap. 

These different parity effects become even more impressive when fixing the 
Zeeman energy and varying the magnetic flux. In the trivial phase, the excitation 
energies for even parity are of order 2 A e g with small oscillations of period $o /2 and 
amplitude d 2 /A c g. For odd parity, the QP gap is absent and the excitation energies 
reflect the single-particle level spacing. In contrast, as shown in figure [ijb) we find 
large $o periodic oscillations of amplitude 2A c ff in the nontrivial phase. Here, the 
excitation energies for &/$o G (—1/4, 1/4) are determined by the superconducting gap 
2A c ff due to the pairwise creation of Bogoliubov QPs i.e. the breaking of Cooper pairs, 
while they display the single-particle level spacing d 2 /A e g for $/$o € (1/4, 3/4). The 
excitation spectrum for odd parity is equivalent to that for even parity but shifted by 
$o/2, as follows from the discussion above. 

This characteristic $0 flux periodicity of the excitation spectrum in the nontrivial 
superconducting phase is directly related to the 4ir periodicity of the Josephson current 
between two topological SCs []I3[33J|35] which has been recently discovered in InSb/Nb 
nanowire junctions [26 . In first order in the tunneling matrix element t, the Josephson 
Hamiltonian between two one-dimensional topological SCs is given by 

Hj{Act>)=Vtcos(^-), (15) 

where V, which has eigenvalues ±1, describes the parity operator of the neutral 
fermion state shared between the two topological SCs and A(j) the superconducting 
phase difference. For a fixed parity of the neutral fermion, the Josephson energy 
is 4-7T periodic in the phase difference. If the Josephson junction is inserted into a 
ring structure, the magnetic flux threading the ring yields a superconducting phase 
difference A<fi = 47r<I>/<f > o, and the An phase periodicity is equivalent to a <&o flux 
periodicity. If the nanowire is coupled to a reservoir, varying A</> by 2tt ~ $0 1 2 will 
change the occupancy (V+l)/2 of the neutral fermion state and hence the parity of the 
ground state. This is in agreement with our result that the parity of the ground-state 
wave function changes when varying the flux through the ring by $0/2. 

3.1. Dependence on details of the geometry 

In this section, we study how details of the geometric realization of the ring topology 
affect the excitation spectrum. For this purpose, we compare the spectra for a ring, a 
square, and for a model with periodic BC as sketched in figure [TJa). In Ref. [23], the 
authors showed that the low-energy physics of the ring-shaped nanowire is equivalent 
to that of a strip of width L± and length L — 2nR 3> L± with periodic BC along the 
^-direction and with vector potential A = (<J> — &o/2)x/L in the Landau gauge. Here, 
the first term in the bracket describes the magnetic flux penetrating the ring. The 
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Figure 1. (Color online) (a) Sketch of three different realizations of the S 1 
topology of the circle: a ring, a square, and a wire with periodic boundary 
condition indicated by dots, (b) Lowest excitation energies for Ez = 2A as 
function of magnetic flux, and (c) lowest excitation energies for $ = (mod <t>o) 
as function of magnetic field. The parameters used in the calculation are: even 
parity, N = 44, L = 2ttR = 3 /im, and a = 2 nm. In (c), the magnetic flux for 
the system with periodic BC is <E> = — $o/2 (mod 2). The lines for panel (b) are 
defined in panel (c). 



second term originates from the conservation of the total angular momentum of the 
electrons in the ring-structure. More specifically, this conservation yields a spin-orbit 
coupling between states with equal total angular momentum (pR, t) and (pR + H, I) 
and thus effectively shifts orbital angular momenta pi? by ±h/2 [49 j. This shift can be 
understood by identifying it with a 2n spin rotation of an electron encircling the ring 
which is equivalent to a Berry phase factor —1 and an effective shift of the magnetic 
flux by -$ /2. 

Experimentally, the fabrication of ring structures with radii of several hundred 
nanometers is challenging and the approximation of the ring by a triangle or a rectangle 
is likely. For this purpose, we consider a square-like structure made out of four 
nanowires, however, our results are qualitatively also valid for a triangular structure 
which consists of three wires. 

In figure [ljb),(c), we compare the fixed particle number excitation spectrum for 
the three different realizations of the lattice on which the Hamiltonian is defined; the 
ring of radius R, a quadratic approximation of the ring with edge lengths ttR/2, and a 
straight nanowire of length L = 2irR with periodic BC. In our numerics, we model the 
ring-shaped nanowire by a one-dimensional tight-binding Hamiltonian Eq. Q with 
spin-orbit coupling perpendicular to the nanowire, and thus rotating in the x-y plane. 
As function of the discretized azimuthal angle tfi — 2ni/n with n lattice sites, the 
spin-orbit direction is then given by 

CTSO,ring(«) = sin(ipi)a x + cos(^)cr y . (16) 

Similarly, we model the square by abrupt changes in the spin-orbit direction at the 
position of the corners, 

a x , for < i < f 
, forf<i<! 

<7SO, S quarcW - \ for | < i < % ( ?j 

-a y , for 5p < i < n 
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and the straight nanowire with periodic BC by a constant spin-orbit direction, 

CSO : pcriodic(«) = <? V ■ (18) 

The spectra for these three models are qualitatively very similar and show a $o 
flux period. As expected, we find that the model with periodic BC yields the same 
spectrum as the ring model but with vector potential shifted by — & x/2L due to the 
Berry phase factor — f. The spectrum for the square model shows small deviations 
from that for the ring geometry since here both translational and rotational symmetry 
are broken. As a consequence, QP states with reduced excitation energy exist, which 
are predominately localized near the corners of the square. Because of the broken 
symmetries, we also find that the spectrum as function of flux is shifted by ~ $0/8 
compared to the ring model. 

We see that our main results are robust against the details of the geometric 
realization and rely on the existence of a void such that the topology of the nanowire 
is homotopically equivalent to an annulus. All these results underline the general 
arguments in the introduction, connecting ground-state degeneracies on the torus to 
parity and flux periodicities of excitations. 



3.2. Electrostatic disorder 

On one hand, disorder is known to often have drastic influence on the electronic 
properties of low-dimensional systems. On the other hand, superconducting pairing 
correlations in s-wave SCs are protected against time-reversal invariant impurity 
scattering by Anderson's theorem [50]. This motivates us to address the question of 
how robust the ground-state degeneracies in the torus topology are against potential 
electrostatic disorder. In the following, we discuss the effect of disorder on the 
excitation energies in the regime where the effective gap is larger than the single- 
particle level spacing A c ff > d, i.e. for 

E Z < *£L. (19) 

We model electrostatic disorder by a locally varying impurity potential V r with 
vanishing mean value and Gaussian white noise correlator (V T V T ') = j8 rjr i/a. We 
here consider the regime of disorder strengths 7 < Ae so Z so , since strong disorder 
7 3> Ae so l so breaks the nanowire into topological and nontopological domain walls 
and thereby destroys the excitation gap |37j . 

In figure [2] we display the excitation spectrum as a function of Zeeman energy 
for different combinations of parity and magnetic flux. We find that the effect of 
electrostatic disorder is very different in the topologically trivial and the nontrivial 
phase. In the trivial phase, the QP excitation gap is remarkably robust against 
disorder which is characteristic for s-wave superconductivity 50J . In contrast, we 
find a significant reduction of the excitation gap due to disorder in the nontrivial 
phase. Since the topological phase for large Zeeman energies Ez ~> A, \i can be 
mapped onto a spinless p-w&ve SC [30], this reduction is in full agreement with the 
effect of disorder on the excitation gap in spinless p-wave SCs [HTJ[52). Furthermore, 
we find that the reduction is very efficient near the topological phase transition since 
there already weak disorder breaks the nanowire into domains of different chemical 
potential and thereby shifts parts of the wire through the topological phase transition 
which reduces the excitation gap locally. Away from the topological phase transition, 
the reduction of the excitation gap is weaker because the existence of partially trivial 
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Figure 2. (Color online) Lowest excitation energy for a ring-shaped nanowirc 
with N = 44, L = 3 fim, and a = 5 nm as function of Zeeman energy and for 
different disorder strengths. The curves represent the average over 50 random 
disorder configurations. The lines are defined in panel (c). 



domains due to disorder becomes unlikely. Furthermore, we find that disorder shifts 
the topological phase transition towards larger values of the Zeeman energy [5H [S3] ■ 
As before, we argue that this shift originates from local topological phase transitions 
at E% = A 2 + (/i + V r ) 2 which are shifted towards larger Zeeman energies due to 
disorder. 

Since the parity and flux dependence of excitation energies reflect the presence 
or absence of nontrivial topological order, our findings for the nonlinear Coulomb 
blockade transport are robust against electrostatic disorder and other perturbations 
as long as the topological order is not destroyed by the formation of domain walls. 

3.3. nonsuperconducting segments 

In this section, we consider the situation that the proximity induced superconducting 
order parameter is spatially dependent. Experimentally this might appear due to the 
roughness of the nanowire/SC interface or if the nanowire is not completely covered 
with the s-wave SC. As sketched in figure ^ a), we describe this spatial dependence of 
the superconducting pairing amplitude by a step function such that A r = for < 
x < A « L and A r = A elsewhere. In figure^c), we display the excitation energies for 
different lengths A of the nonsuperconducting segment. We find a significant reduction 
of the excitation gap in the trivial phase while the excitation energies in the nontrivial 
phase are only weakly reduced even for A ~ l so . We argue that the origin of the 
robustness of the excitation gap in the topological phase is the small effective gap 
A e ff < A and the enhanced Fermi velocity vp(Ez > A) ~ 2vp(Ez = 0) due to the 
occupation of a single spinless band. Hence, the superconducting coherence length 
£ = vp/A C H in the topologically nontrivial phase is significantly enhanced as compared 
to the trivial phase. Thus, in the nontrivial phase superconducting pairing correlations 
are more efficiently induced in the nonsuperconducting part of the wire which here 
shows up as the robustness of the superconducting gap against the existence of a quite 
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Figure 3. (Color online) Sketch of the nanowire with (a) one and (b) ten 
nonsuperconducting segments of length A, i.e. segments which are not covered 
with superconducting material, (c), (d) Lowest excitation energy for a ring-shaped 
nanowire of length L = 3 fim with even parity, N = 44, and <t> = for the case 
with (c) one and (d) ten nonsuperconducting segments. 



long normal segment. 

From the robustness of the excitation gap, we conclude that it is not necessary 
for our proposed setup that the nanowire is completely covered with the s-wave SC. 
In particular, we propose that it is sufficient to place superconducting grains on the 
nanowire in order to significantly increase the charging energy and to reduce Cooper 
pair cotunneling through the SC [see section [5]. We now assume that the nanowire 
contains ten nonsuperconducting segments of length A uniformly distributed over the 
nanowire as sketched in figure [3jb). In figure [3jcL) , we study the excitation spectrum 
for different characteristic lengths A. While the excitation gap for A = 3l so /2 ~ 150 
nm in the trivial phase is completely absent, we find that the excitation energies in the 
nontrivial phase are only reduced by 30% as compared to the situation where A ^ 
everywhere. 

4. Multi-band Hamiltonian 

In this section we make a departure from the case of strictly one-dimensional nanowires 
and consider the experimentally realistic situation of quasi one-dimensional nanowires 
of finite thickness with a <C Lj_ < £. In our numerical analysis, we model the ring 
shaped nanowire by a strip with periodic boundary conditions along the x-direction 
and with hard wall boundary conditions along the y-direction. The magnetic flux $ is 
incorporated through the modified vector potential A = ($ — §q/2)xJL as discussed 
above in section [XT] 

To ensure that the induced superconducting phase remains quasi one-dimensional 
and the nanowire exhibits a substantial gap, we demand that the width does not exceed 
the superconducting coherence length £ = u_p/A e g [MJ [55J ESI [53 EE] • The spatial 
extension in the y-dircction gives rise to the existence of additional transverse modes 
and thus subbands which might be partially occupied depending on the chemical 
potential. In figure |4^a), we display the Bogoliubov QP spectrum for $ = — $o/2 as 
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Figure 4. (Color online) (a) Lowest energy of the Bogoliubov QP spectrum 



phase diagram inferred from the topological number Q as function of Zeeman 
energy and mean electron number. The green dashed lines in (a) represent the 
topological phase transition for the various subbands in the limit of uncoupled 
subbands a — > 0, whereas the color scale encodes the excitation energy for coupled 
subbands. The parameters used in the calculation are: L = 3 fim, L± = 100 nm, 
a = 5 nm, and 3> = — <t>o/2. 



function of Zeeman energy and chemical potential. For fi < (irti) 2 /2mL 2 L only one 
subband is partially occupied and the excitation spectrum is equivalent to the one 
discussed in section [3j With increasing chemical potential higher subbands are filled 
up consecutively and similarly to the single-band case, the higher subbands can be 
either topologically trivial or nontrivial depending on the chemical potential and the 
Zeeman energy. 

In figure [4| a), the topologically nontrivial phase shows up as islands which are 
enclosed by lines of vanishing excitation energies, i.e. by topological phase transitions. 
Assuming that the subbands are uncoupled, we find the topological phase whenever 
the chemical potential lies well within one of the spin gaps at zero momentum and 
when the Zeeman energy satisfies the relation 

E% > A 2 + (fi-e n ) 2 , (20) 

where e„ = (hnir) 2 /2mL 2 L denotes the kinetic energy of subband n. However, 
the transverse spin-orbit term ac\ (J a^ ra ,c r+ s y a i couples the subbands and thereby 
modifies the lines where the topological phase transitions occur. These modifications 
are similar to avoided crossings with energy splitting <5/i rs 2akp lV between the lines of 
topological phase transitions and thus reduce the size of the topologically nontrivial 
islands. With increasing chemical potential, the spin-orbit energy in the transverse 
direction increases and thus the energy splitting due to the avoided crossing increases 

In figure[4](b), we display the topological number Q as function of Zeeman splitting 
and mean electron number N. As before, the topological number is +1 in the trivial 
and —1 in the nontrivial phase and thus we conclude that the parity of the grand 
canonical ground states in both phases is different with even parity in the trivial 
and odd parity in the nontrivial phase. Similarly to figure ^ a), we find islands of 
topologically nontrivial phase which are enclosed by the trivial phase. We propose 
that the fixed mean particle number excitation energies can be used as a tool to 
investigate the topological phase diagram. In |Appendix A we show that our results 
for the single-band nanowire can directly be applied to the multi-subband nanowire. 
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Figure 5. (Color online) Lowest excitation energy for fixed mean particle number 
as function of magnetic flux <t> (mod <t>o) which pierces the nanowire relative to its 
central line, and as function of additional magnetic flux <5<t> = BLL±/2 through 
the nanowire itself due to the finite thickness of the nanowire. The parameters 
used in the calculation are: odd parity, N = 43, Ez = 2A, L = 2irR = 3 fim, 
Lj_ = 70 nm, and a = 10 nm. For even parity, we find the same spectrum with <t> 
shifted by <t> /2. 



Thus, the topologically trivial phase is characterized by an excitation spectrum with 
gap 2A c ff for even parity while the excitations are determined by the single-particle 
level spacing for odd parity. When varying the magnetic flux, the spectra for both 
even and odd parity show small $o/2 periodic oscillations as expected for trivial SCs. 
In contrast, the situation is different in the topologically nontrivial phase where the 
excitation spectrum qualitatively depends on both magnetic flux and electron parity. 
We here find a characteristic 4>o flux period similar to the situation for the single- 
band model in section [3] We conclude that the excitation spectrum for fixed mean 
particle number, which can be observed in nonlinear Coulomb blockade transport, is 
an unbiased tool to map out the topological phase diagram shown in figure Bib). 




Due to the finite width Lj_ of the nanowire, the area of the nanowire itself is 
penetrated by magnetic flux and thus the magnetic flux through the ring-shaped 
nanowire is not well-defined. However, the magnetic flux can be decomposed into a 
mean value for the middle of the wire and deviations due to the finite thickness 



For nanowires with radius R = L/2n ~ 0.5 /J,m and magnetic field strengths B < 1 T, 
we find (5<I>(L^/2) > $0- In figure [5j we display the fixed electron number excitation 
spectrum as function of mean magnetic flux $ (mod $ ) an d additional flux <5$(Lj^/2). 
We find that the flux periodicity of the excitation spectrum is not changed, however, 
the excitation spectrum itself is shifted due to <5$. We can therefore conclude that 
additional magnetic flux due to the finite width of ring-shaped nanowires with large 
radii is unproblematic for the study of the flux periodicity of the excitation energies. 

5. Sequential and Cotunneling of Cooper pairs 

In general, there are competing transport channels through the hybrid ring-shaped 
nanowire. Our focus is on nonlinear Coulomb blockade transport due to sequential 
single-electron tunneling through the SM. The most important competing channels 
are sequential and cotunneling of Cooper pairs. For that purpose, we estimate the 




5§(y) = BLy for 




(21) 
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magnitude of the current due to Cooper pair processes for a general superconducting 
island weakly tunnel-coupled to two leads and study under which parameter conditions 
they become important. We assume that the SC has a large number of transverse 
channels N± and charging energy E c . The tunneling between lead i and the SC is 
described by 

H T ,i = { < k, q aL CT ( M q CT 7 qCT + *V7- q - CT ) + h.c.y (22) 

kqfT 

where ik q & r e the tunnel matrix elements, are the fermion operators in lead i 
with energy eik, and 7 q(T are QP operators for the SC with energy Eq = w£2 -). /\| 



The Uqa and i> qcr are the BCS coherence factors with magnitudes \/(l ± £q/Eq)/2. 

Tunneling of Cooper pairs between lead i and the SC is described by an effective 
Hamiltonian which can be derived in second order perturbation theory in Ht,%- In the 
first step, one electron with momentum and spin a is transferred from an initial 
state into an intermediate excited state with momentum q of the superconducting 
island. In the second step, another electron with momentum k2 and spin — a tunnels 
into the partner state of the first electron — q such that both electrons form a Cooper 
pair. Hence, the final state contains an extra Cooper pair in the SC and two QP 
excitations in the lead. Similarly, we find the reverse process by splitting a Cooper 
pair followed by two consecutive electron tunneling events [59J • This yields the effective 
tunneling Hamiltonian 

H C p, = (BCS|g T . a . 1 „ ffr, t |BCS), (23) 
irj - H 

where we traced out the QP operators via the BCS ground state |BCS). We find 

Hcpa = E {^k 1 ,k 2 aikitaik 2 4 + ^*k 1 ,k 2 ct lk 2 j. a lk 1 t} ( 24 ) 

kik 2 

with the effective tunneling matrix elements 

1 



+ F t F - f "}• (2r»i 

In the following, we consider the Andreev current through a normal- 
superconducting-normal structure with symmetric barriers and bias voltage < V < 
Asc/e. Assuming that the voltage between the left (right) lead and the SC is ±V/2, 
we calculate the rate for the Andreev reflection process using Fermi's golden rule. The 
current for the scattering of two electrons from the left metallic lead into the SC reads 

I A (u>) = 2e— J2 \AL kl ,k 2 \ 2 f(e Lkl )f(eL k M£L kl +e Lk2 +cj) (26) 

kik 2 

with u = eV — 4E C and Fermi functions /. In Ref. [59] it has been shown that the 
Andreev conductance Ga — Ia/V f° r sequential Cooper pair tunneling can be written 
as Ga(w) = (e 2 /h)G 2 (oj)/N±, where G is the dimensionless normal state conductance 
and A^_l the number of transverse channels through the superconducting region in its 
normal state. Due to the charge 2e of Cooper pairs, sequential tunneling of Cooper 
pairs is not resonant for eV/2 < E c — A c g and can be neglected. In the expression for 
the current, this suppression shows up as a shifted chemical potential uj — eV — 4E C . 
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Similarly, we calculate the current for Cooper pair cotunneling from the left lead 
to the right lead via the superconducting island by calculating the scattering rate in 
second order perturbation theory in Hcp.i- We find 

I Acot (V) = 2e 87T V- \A Lkl M 2 \ 2 \A Rka 



H u Ct u ( £ ik! + £Lk 2 +eV - AE C ) 2 

kik 2 k 3 k 4 



/(eLk 1 )/(eLk 2 )/(-ei?k 3 )/(-eflk 4 ) 

^(eiki + £ik 2 ^ £i?.k 3 - £-Rk 4 + 2eV). (27) 

Building on the result for the sequential Cooper pair tunneling and assuming eV < -Ec, 
we find that the Andreev cotunneling current can be expressed as 

(28) 

Et 



In the expression for Eq. (28), the Andreev conductance Ga(cV) is not suppressed by 
the Coulomb energy since the charge on the superconducting island after the tunneling 
events is the same as the initial charge. 

In contrast, we find for sequential electron tunneling a current I seq — (e/h)T 
where Y is the tunneling rate between the lead and the SM. For characteristic bias 
voltages smaller or equal to E c /e, we compare the currents due to the sequential 
tunneling of electrons and the Andreev cotunneling of Cooper pairs. With Eq. (28) 
and the expression for the Andreev conductance, we find 

J^- « (29) 

We now make the conservative assumption T w d/10 and d ss E c /10, where d is the 
mean level spacing in the SM, and demand that single particle sequential tunneling 
be larger than Cooper pair cotunneling. In this way, we obtain the condition that 
G < \/N_i_/3, i.e. the dimensionless conductance of the junction between lead 
and the SC in its normal state has to be smaller than one third of the square 
root of the number of transverse channels. For a metal of diameter 10 nm and 
with Fermi wavelength 0.3 nm, the number of transverse channels is approximately 
(diameter/wavelength) 2 = 1000, and thus the dimensionless normal state conductance 
needs to satisfy G < 10, which is realistic for metallic quantum dots with current state 
technology. 

One way to realize the condition G < yN j_ /3 experimentally is to not cover 
the nanowire with superconducting material in the vicinity of the electrodes. This 
significantly reduces the conductance between the SC and the electrodes. One can 
even imagine that an extreme limit could be realized, in which all electrons entering 
the hybrid system have to do so via the SM in the vicinity of the electrodes. One might 
argue that as a consequence of removing the SC near the electrodes, the proximity 
induced pairing amplitude in this region will be reduced as well. However, when the 
region not covered with superconducting material is considerably smaller than the 
coherence length in the SM (of the order of 100 nm as shown in Sec. 3.3), this effect 
will be small. In principle, one could go even further and only deposit superconducting 
nanograms on top of the nanowire instead of adding a fully connected SC, and in this 
way eliminate the influence of Andreev cotunneling almost completely. 

In order to fully suppress cotunneling of Cooper pairs through the SC, we propose 
to use ferromagnetic leads with the polarization in magnetic field direction. While 
ferromagnetic leads fully suppress Andreev processes and thus cotunneling of Cooper 
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pairs in conventional s-wave SCs, they do not significantly affect the current due to 
sequential tunneling of electrons. 

6. Summary 

In conclusion, we have proposed a Coulomb blockade transport experiment to 
investigate the topological order of semiconductor-superconductor hybrid nanorings, 
and have shown that characteristic parity and flux periodicity effects in the excitation 
spectrum reflect the distinct ground-state degeneracies of trivial and nontrivial 
superconducting phases on manifolds with nonzero genus. In particular, the excitation 
spectrum for fixed mean particle number provides clear signatures of the h/e flux 
periodicity in the nontrivial phase and the topological phase transition. All these 
findings are robust against geometry details of the realization of the ring structure 
and rely on the existence of a hole such that the system is homotopically equivalent 
to a circle. 

We have shown that the spectroscopic gap in the nontrivial phase is robust against 
moderate electrostatic disorder. Furthermore, the nontrivial phase is characterized 
by a large superconducting coherence length which allows to deposit superconducting 
nanograins on top of the nanowire instead of adding a fully connected superconductor, 
and in this way reduces the Andreev cotunneling and enhance the charging energy. 
Using a T-matrix formalism, we have estimated the magnitude of Andreev cotunneling 
and have derived a criterion for the maximum number of parallel conduction channels 
through the proximity coupled s-wave superconductor which ensures that single 
particle transport dominates over cotunneling of Cooper pairs. 

Finally, we studied multi-subband nanowires and we have shown that nonlinear 
Coulomb blockade transport can be used as a tool to map out the topological phase 
diagram. 
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Appendix A. Multi-band phase diagram 

In this appendix, we present the lowest excitation energy En — E 9 N S for the multi-band 
SM hybrid nanowire for several combinations of magnetic flux and parity as function 
of Zeeman energy and mean electron number. As shown in figure [AT] and explained in 
section [4j both the chemical potential pL and the Zeeman energy can be used to tune 
the nanowire through the topological phase transitions. 

Our findings for the trivial phase (i.e. the dark region in figure |Al[ a)) 
are characteristic for s-wave superconductivity in superconducting grains without 
excitation gap for odd parity [figures |AT| a) and (c)] and with energy gap 2A e g for even 
parity [figures |AlJb) and (d)]. These excitation energies do not change qualitatively 
when changing the magnetic flux and show small < I ) o/2 periodic oscillations of order 
d 2 /A eS < A off . 

In the topologically nontrivial phase (i.e. the bright region in figure [Alja)) the 
parity effect is very different. Here, the excitation energies depend on both electron 
parity and magnetic flux. In Figs. |Al[ a) and (d) we find an excitation gap 2A c ff 
since two Bogoliubov QP excitations are required and thus a Cooper pair needs to be 
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Figure Al. (Color online) Lowest excitation energy in the fermionic excitation 
spectrum with fixed mean electron number as function of magnetic field and mean 
electron number for several combinations of magnetic flux and electron parity; 
L = 3 fim, L± = 100 nm, and a = 5 nm. Note the different color scale in (c) 
where all excitation energies are determined by the single-particle level spacing d. 



broken. In contrast, the excitation energies in Figs. |Al[ b) and Al ]c) are determined by 
the single-particle level spacing since always one unpaired particle is located near the 
Fermi surface. As shown in Fig. |Al[ d), the characteristic signature of the topological 
phase transition is the closing and reopening of the excitation gap. When studying 
the flux dependence of the excitation energies in the nontrivial phase, we find large 
oscillations with period $o and amplitude 2A c ff. For even parity, the excitation 
energies for $/<I>o G ( — 1/4, 1/4) are determined by the effective gap 2A e g while they 
are determined by the single-particle level spacing d 2 / A e g for Q/Qq E (1/4, 3/4). For 
odd parity, we qualitatively find the same spectrum but shifted by $o/2, as follows 
from the earlier discussion. 

Thus, the excitation spectrum for fixed electron number directly reflect the 
topological phase diagram shown in figure |4|b) . The proposed nonlinear Coulomb 
blockade transport experiment can therefore be used as a tool to clearly determine the 
topological order of the hybrid system by measuring the fermionic excitation spectrum. 
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